Perfect state transfer over distance-regular 

spin networks 

M. A. Jafarizadeh a ' 6 ' c *and R. Sufiani ' 6 t 

a Department of Theoretical Physics and Astrophysics, University of Tabriz, Tabriz 51664, Iran. 
^Institute for Studies in Theoretical Physics and Mathematics, Tehran 19395-1795, Iran. 
c Research Institute for Fundamental Sciences, Tabriz 51664, Iran. 

February 1, 2008 



* E-mail: j af arizadeh® t abrizu . ac . ir 
tE-mail:sofiani@tabrizu.ac.ir 



1 



Perfect state transfer 



2 



Abstract 

By considering distance-regular graphs as spin networks, first we introduce some 
particular spin Hamiltonians which are extended version of those of Refs.[l, 2]. Then, 
by using spectral analysis techniques and algebraic combinatoric structure of distance- 
regular graphs such as stratification introduced in [4, 5] and Bose-Mesner algebra, we 
give a method for finding a set of coupling constants in the Hamiltonians so that a par- 
ticular state initially encoded on one site of a network will evolve freely to the opposite 
site without any dynamical controls, i.e., we show that how to derive the parameters of 
the system so that perfect state transfer (PST) can be achieved. As examples, the cycle 
networks with even number of vertices and <i-dimensional hypercube networks are consid- 
ered in details and the method is applied for some important distance-regular networks 
in appendix. 

Keywords: Perfect state transfer, Spin networks, Association scheme, 
Stratification, Distance-regular network 
PACs Index: 01.55.+b, 02.10.Yn 
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The transfer of a quantum state from one part of a physical unit, e.g., a qubit, to another part 
is a crucial ingredient for many quantum information processing protocols [6]. Currently, there 
are several ways of moving data around in a quantum computer. While some methods transfer 
quantum states by moving them down a linear array of qubits, there are others which exploit 
the quantum property of entanglement for teleporting quantum states between distant qubits 
[7, 8]. In a quantum-communication scenario, the transfer of quantum states from one location 
A to another location B, is rather explicit, since the goal is the communication between distant 
parties A and B ( e.g., by means of photon transmission). Equally, in the interior of quantum 
computers good communication between different parts of the system is essential. The need is 
thus to transfer quantum states and generate entanglement between different regions contained 
within the system. There are various physical systems that can serve as quantum channels, 
one of them being a quantum spin system. This can be generally defined as a collection of 
interacting qubits (spin-1/2 particles) on a graph, whose dynamics is governed by a suitable 
Hamiltonian, e.g., the Heisenberg or XY Hamiltonian. 

Quantum communication over short distances through a spin chain, in which adjacent 
qubits are coupled by equal strength has been studied in detail, and an expression for the 
fidelity of quantum state transfer has been obtained [9, 10]. Similarly, in Ref. [11], near perfect 
state transfer was achieved for uniform couplings provided a spatially varying magnetic field 
was introduced. The propagation of quantum information in rings has been also investigated 
in [12]. In our work we focus on the situation in which state transference is perfect, i.e., the 
fidelity is unity, and in which we can design spin networks such that this can be achieved 
over arbitrarily long distances. We will also consider the case in which no external control 
is required during the state transference, i.e., we consider the case in which we have, after 
manufacturing the network, no further control over its dynamics. In general this will lead us to 
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think about more complicated spin networks than the linear chain or chains with preengineered 
nearest-neighbor interaction strengths. We provide two alternative methods for understanding 
how perfect state transfer is achieved with preengineered couplings. This paper expands and 
extends the work done in [1, 2]. We will consider distance-regular graphs as spin networks 
in the sense that with each vertex of a distance-regular graph a qubit or a spin is associated 
(although qubits represent generic two state systems, for convenience of exposition we will 
use the term spin as it provides a simple physical picture of the network). Then, duo to 
the fact that distance-regular graphs are underlying graphs of association schemes (see for 
example [13, 20]), we use their algebraic properties in order to find suitable coupling constants 
in some particular spin Hamiltonians so that perfect transference of a quantum state between 
antipodes of the networks can be achieved. More clearly, for a given distance-regular network 
first we stratify the network with respect to an arbitrary chosen vertex of the network called 
reference vertex (for details about stratification of graphs, see [4, 5, 20]). Then, we consider 
coupling constants so that vertices belonging to the same stratum with respect to the reference 
vertex possess the same coupling strength with the reference vertex whereas vertices belonging 
to distinct strata possess different coupling strengths. Then we give a method for finding a 
suitable set of coupling constants so that PST over antipodes of the networks be possible. As 
examples we will consider the cycle networks with even number of vertices and d- dimensional 
hypercube networks in details and some important distance-regular networks in an appendix. 

The organization of the paper is as follows: In section 2, we review some preliminary 
facts about association schemes, stratification, distance-regular graphs and spectral analysis 
techniques. Section 3 is devoted to perfect state transfer (PST) over antipodes of distance- 
regular networks, where a method for finding suitable coupling constants in particular spin 
Hamiltonians so that PST be possible, is given. The paper is ended with a brief conclusion 
and an appendix. 
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In this section we will first review some preliminary facts about distance-regular graphs, cor- 
responding stratification and spectral distribution associated with graphs. 

2.1 Association schemes 

First we recall the definition of association schemes. For further information on association 
schemes, the reader is referred to Ref. [13]. 

Definition. An association scheme with d associate classes on a finite set V is a set of matrices 
A , A 1 , ...,A d in R VxV , all of whose entries are equal to or 1, such that 

(i) A = I V ; 

(ii) Ai is symmetric for i — 1, d] 

(iii) for all i,j in {0, 1, d}, the product AiAj is a linear combination of A , A x , A d ; 

(iv) none of the Ai is equal to O v , and Yn=o Ai — J v , where v := \ V\ and J v is an v x v all one 
matrix. 

It should be noticed that, since Ai is symmetric with entries in {0, 1}, the diagonal entries of 
Af are the row-sums of A^ Condition (iii) implies that Af has a constant element, say on 
its diagonal. Therefore every row and every column of Ai contains «j entries equal to 1. Hence 
Ai J v = J v Ai = KiJ v . Moreover, A Ai = AiA = Ai. 
From condition (iii), one can write 

^ = X>*A, (2-1) 

which implies that the adjacency matrices A , A±, A d form a basis for a commutative algebra 
A known as Bose-Mesner algebra associated with the association scheme. This algebra has a 
second basis E ,...,E d such that, EiEj = 5ijEi and X^Lo &i = I with E = l/vJ v [13]. The 
matrices Ei for < i < d are known as primitive idempotents of Y . Furthermore, there are 
matrices P and Q such that the two bases of the Bose-Mesner algebras can be related to each 
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other as follows 



d 



3=0 



1JL 



(2-2) 



3=0 



Then, clearly we have 



PQ = QP = vl. 



(2-3) 



It also follows that 



AjEi = 



P E 



(2-4) 



which indicates that Pjj is the 2-th eigenvalue of A,- and that the columns of E i are corre- 
sponding eigenvectors. Also, m 8 := trEi = v(a\E i \a) = Q i0 (where, we have used the fact that 
(a\Ei\a) is independent of the choice of a E V, see Eq.(2-3)) is the rank of the idempotent Ei 
which gives the multiplicity of the eigenvalue P^ of Aj (provided that P^ ^ Pkj for k ^ i). 
Clearly, we have m = 1 and Ya=q frii — v since Yh=o Ei — I. 

The underlying network of an association scheme T = (V, E) is an undirected connected 
network with adjacency matrix A = A\. Obviously replacing A\ with one of the other adja- 
cency matrices A iy i ^ 0, 1 will also gives us an underlying network V = {V,E') (not necessarily 
a connected network) with the same set of vertices but a new set of edges. 

As we will see in subsection 2.3, in the case of distance-regular networks, the adjacency 
matrices Aj are polynomials of the adjacency matrix A = A 1 , i.e., Aj = Pj(A), where Pj is a 
polynomial of degree j, then the eigenvalues P^ in (2-4) are polynomials of eigenvalues Pa = \ 
(eigenvalues of the adjacency matrix A). This indicates that in distance-regular graphs, the 
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matrix P* is a polynomial transformation [14] as 

1 1 ... 1 

Pi(Xo) Pi(Xi) ... Pi(X d ) 

P 2 (A ) P 2 (Ai) ... P 2 (A d ) 



^ P d (A ) P d (A0 . . . P d (X d ) ) 



(2-5) 



or P Jl = P(A,). 



2.2 Stratifications 

For a given vertex a 6 V, let T^a) := {(3 G V : (a,/3) G P«} denotes the set of all vertices 
having the relation Pj with a. Then, the vertex set V can be written as disjoint union of 
for i — 0, 1, 2, d, i.e., 

\/=[Jr,(a). (2-6) 

i=0 

We fix a point o G V as an origin of the underlying graph of an association scheme, called 
reference vertex. Then, the relation (2-6) stratifies the graph into a disjoint union of associate 
classes Tj(o) (called the i-th stratum with respect to o). Let 1 2 {V) denote the Hilbert space 
of C-valued square-summable functions on V. With each associate class Tj(o) we associate a 
unit vector in l 2 (V) defined by 

\k) = A= E l«>> (2-7) 

V K i aeTi(o) 

where, \ct) denotes the eigenket of a-th vertex at the associate class Tj(o) and k, l = |Fj(o)| 
is called the i-th valency of the graph. Now, let A4 be the adjacency matrix of the graph 
T = (V, P). Then, for the reference state \<j> ) (|0 O ) = |o), with o G V as reference vertex), one 
can write 

A|0o>= E (2-8) 

/3 G r,(o) 
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Then, by using (2-7) and (2-8), we obtain 

AM = V^M- (2-9) 

One should notice that, in underlying networks of association schemes, stratification is reference 
state independent, namely one can choose any arbitrary vertex as a reference state. 

2.3 Distance-regular networks and spectral techniques 

Distance-regular graphs are underlying graphs of so called P-polynomial association schemes 
[13], where the adjacency matrices A^ are defined based on shortest path distance. More 
clearly, if distance between nodes a, (3 G V denoted by d(a, (3) be the length of the shortest walk 
connecting a and (3 (recall that a finite sequence a , a±, a n G V is called a walk of length n if 
~ «fc for all k — 1, 2, n, where ctk-i ~ ctk means that a^-i is adjacent with then the 
adjacency matrices Ai for i — 0, 1, d in distance-regular graphs are defined as: (A)a,/3 = 1 
if and only if d(a,f3) = % and (A)a,/3 — otherwise, where <i :=max{<9(o;, /5) : a,/? 6 V} is 
diameter of the graph. 

For distance-regular graphs, the non-zero intersection numbers are given by 

Oi=p\i, h=pl +ll , Ci=p]_ hl , (2-10) 

respectively. The intersection numbers (2-10) and the valencies Ki with K\ = k{= deg(a), for 
each vertex a) satisfy the following obvious conditions 

a i ~\~ bi + Cj = /t, Ki—ibi—i = /tjCj, i = 1, <i, 

/to = ci = 1, &() = «;! = /c, (c = b d = 0). (2-11) 

Thus all parameters of a distance-regular graph can be obtained from its intersection array 
{b , bd-i, Ci, Cd}- Then, it can be shown that the following recursion relations are satisfied 

A\A,i = bi-xA^x + ciiAi + c i+1 A i+1 , i = 1, 2, d - 1, 
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A x A d = & d _i4j-i + (« - c d )A d . (2-12) 
The recursion relations (2-12) imply that 

A i = P i (A), i = 0,l,...,d. (2-13) 

Then, one can easily obtain the following three term recursion relations for the unit vectors 
\</>i), % = 0, 1, ...,d 

A\<j>i) = A+il^+i) + + A|0i-i>, (2-14) 

where, the coefficients «j and $ are defined as 

a = 0, a k = a k = K — b k — c k , uj k = 01 = & fe _ic fe , k = \,...,d, (2-15) 

(see Ref. [5, 21, 22] for more details). 

Now, we recall some preliminary facts about spectral techniques used in the paper, where 
more details have been given in Refs. [5, 21, 22]. 

For any pair (A, \(f>o)) of a matrix A and a vector \<f>o), one can assign a measure \i as follows 

u.(x) = (cf) \E(x)\ ( f) ), (2-16) 

where E(x) = J2i \ui)(ui\ is the operator of projection onto the eigenspace of A corresponding 
to eigenvalue x, i.e., 

A = J xE(x)dx. (2-17) 
Then, for any polynomial P(A) we have 

P(A) = J P(x)E(x)dx, (2-18) 

where for discrete spectrum the above integrals are replaced by summation. Therefore, using 
the relations (2-16) and (2-18), the expectation value of powers of adjacency matrix A over 
reference vector |0 O ) can be written as 

(0 o |A m |0 o ) = / x m fi(dx), m = 0, 1, 2, .... (2-19) 

JR 
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Obviously, the relation (2-19) implies an isomorphism from the Hilbert space of the stratifica- 
tion onto the closed linear span of the orthogonal polynomials with respect to the measure \i. 
More clearly, from orthonormality of the unit vectors \<pi) (with |0 O ) as unit vector assigned to 
reference node) we have 

Sij = (<f>i\<f>j) = -^(0o|A^# o > = / Pl(x)P^x)^dx), (2-20) 

with P((A) := -j=Pi(A) where, we have used the equations (2-9) and (2-13) to write 

A\<M = ^=Pi(A)\<M = P' i (A)\<l Kt ). (2-21) 



Now, by substituting (2-21) in (2-14) and rescaling P' k as Q k = (3i...f3 k P' k) the spectral distribu- 
tion n under question will be characterized by the property of orthonormal polynomials {Q k } 
defined recurrently by 

Qo{x) = 1, Qi(x) = x, 
xQ k (x) = Q k+1 (x) + a k Q k (x) + PlQ^x), k > 1. (2-22) 
where, the coefficients ctj and $ are defined as 

a k = n-b k -c k , u k = Pi = b k _!C k , k = l,...,d, (2-23) 

If such a spectral distribution is unique, the spectral distribution ji is determined by the 
identity 

G ^ Lx-y x _ aQ Qd+i(x) ^i-ii' 



x — an 

z x~ CK3 

where, x t are the roots of polynomial Q d+ i(x). G^(x) is called the Stieltjes/Hilbert transform 
of spectral distribution \x or Stieltjes function and polynomials {Q k ^} are defined recurrently 

as 

Q { o\x) = 1, Q?\x) = x - ai , 
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xQt\x) = QiUx) + a k+l Qi\x) + Pt^QiUx), k > 1, 



(2-25) 



respectively. The coefficients 7; appearing in (2-24) are calculated as 



7i := lim [(x - Xi)G^(x)} 



(2-26) 



Now let G^(z) is known, then the spectral distribution \i can be recovered from G^(z) by 
means of the Stieltjes/Hilbert inversion formula as 



Substituting the right hand side of (2-24) in (2-27), the spectral distribution can be determined 
in terms of Xi, I — 1, 2, ... and Guass quadrature constants 7;, / = 1, 2, ... as 



(for more details see Refs.[16, 17, 18, 19]). 

3 Perfect State Transfer (PST) over antipodes of distance- 
regular networks 

3.1 State Transfer in Quantum Spin Systems 

The PST algorithm was proposed by Christandl et al. [1,2], and it can be implemented in the 
XY chain. The algorithm can transfer an arbitrary quantum state between the two ends of 
the chain in a fixed period time, only using XY interactions. For one-dimensional fermionic 
chains, the model of a system consisting of spinless fermions (or bosons) hopping freely in a 
network of N lattice sites can be mapped to spin chains in which spins are coupled through 
the XY Hamiltonian 




(2-27) 



(2-28) 



H = - 



1 N-l 1 N 



E JA°PU + ° v i°U) + E A >J + !)> 



(3-29) 




Perfect state transfer 



12 



by the Jordan- Wigner transformation, where Jj is the time-independent coupling constant 
between nearest-neighbor sites j and j + 1, and \j represents the strength of the external 
static potential at site j. 

A quantum spin system associated with a simple, connected, finite graph G = (V, E) as a 
spin network is defined by attaching a spin-1/2 particle to each vertex of the graph so that to 
each vertex i e V one can associate a Hilbert space Hi ~ C 2 . The Hilbert space associated 
with G is then given by 

n G = ® iev n l = (c 2 r N , (3-30) 

where N := \V\ denotes the total number of vertices in G. On the other hand, quantum state 
transfer over a network is similar to the quantum random walk problem, where a variety of 
networks are equivalent to one-dimensional chains [1,22]. Therefore, it can be focused on a 
chain of N sites. For j = 1,2, ...,N, let \j) be the state where a single fermion (or boson) 
is at the site j but is in the empty state |0) for all other sites and |0) be the vacuum state 
where all sites are empty. For spin chains, |0) corresponds to the state where all the spins are 
in the spin-down state | |) and \j) corresponds to a spin-up state | f) f° r the jth spin and 
spin-down for all other spins. The Hamiltonian in this single-particle subspace can be written 
in a tridiagonal form, which is real and symmetric: 



Ai 


Jl 


... 





\ 


Ji 


A 2 


J 2 ... 










J2 


A 3 ... 













Jn-i 










Jn-i 


Xn 


/ 



The quantum state transfer protocol involves two steps: initialization and evolution. First, a 
quantum state \iP)a = ci\0)a+P\1)a £ T~La (with cc, f3 e C and |a| 2 + \f3\ 2 = 1) to be transmitted 
is created. The state of the entire spin system after this step is given by 

\iP(t = 0)) = |^0.. .00 B ) = a\0 A 0...00 B ) + (3\1 A 0...00 B ) = a\0) + (3\A), (3-32) 
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with |0) := |OaO...OOb). Then, the network couplings are switched on and the whole system is 
allowed to evolve under U(t) = e~ tHt for a fixed time interval, say t . The final state becomes 

N 

mo)) = a\0)+(3j2fM\j) (3-33) 
j'=i 

where, fjA(to) '■= {j\e~ lHto \A}. Any site B is in a mixed state if |/ab(£o)| < 1, which also 
implies that the state transfer from site A to B is imperfect. In this paper, we will focus only 
on perfect state transfer. This means that we consider the condition 

\fAB(to)\ = 1 for some < t < oo (3-34) 

which can be interpreted as the signature of perfect communication (or perfect state transfer) 
between A and B in time t . The effect of the modulus in (3-34) is that the state at B, after 
transmission, will no longer be but will be of the form 

a|0) + e^/3|l). (3-35) 

The phase factor e 1 ^ is not a problem because is independent of a and (5 and will thus be 
a known quantity for the graph, which we can correct for with an appropriate phase gate (for 
more details see for example [1, 2, 3]). 

The model we will consider is a distance-regular network consisting of iV sites labeled by 
{1, 2, N} and diameter d. Then we stratify the network with respect to a chosen reference 
site, say 1, and assume that the network contains only the output site N in its last stratum 
(i.e., \(pd) = \N)). At time t — 0, the qubit in the first (input) site of the network is prepared 
in the state \ipi n ). We wish to transfer the state to the Nth (output) site of the network with 
unit efficiency after a well-defined period of time. Although our qubits represent generic two 
state systems, for the convenience of exposition we will use the term spin as it provides a 
simple physical picture of the network. The standard basis for an individual qubit is chosen 
to be {|0) = | I), |1) = | t)}, and we shall assume that initially all spins point "down" along 
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a prescribed z axis; i.e., the network is in the state |0) = |0a00...00b). Then, we consider the 
dynamics of the system to be governed by the quantum-mechanical Hamiltonian 

H G = \jlJ m E Ha, (3-36) 

m=0 (i,j)£Rm 

with Hij as 

Hij = Ci ■ Cj, (3-37) 

where, cr, is a vector with familiar Pauli matrices of, erf and erf as its components acting on 
the one-site Hilbert space Hi, and J m is the coupling strength between the reference site 1 and 
all of the sites belonging to the m-th stratum with respect to 1. 

The total spin of a quantum-mechanical system consisting of N elementary spins <?j on a 
one-dimensional lattice or better called chain is given by: 

N 

(3-38) 

1=1 

One can easily see that, the Hamiltonian (3-36) commutes with the total Spin operator (conser- 
vation). That is, since the total z component of the spin given by o z tot = J2iev °t * s conserved, 
i.e., [<7t ot ,Hc] = 0, hence the Hilbert space Tic decomposes into invariant subspaces, each of 
which is a distinct eigenspace of the operator of ot (this property would be important to use its 
symmetry to diagonalize the Hamiltonian in the well known Bethe ansatz approach). 

In order to consider perfect quantum state transfer, we write the hamiltonian (3-36) in 
terms of the adjacency matrices Ai, i = 0, 1, d of the underlying graph in order to use the 
techniques introduced in section 2 such as stratification and spectral distribution associated 
with the graph. To do so, we recall that the kets \ii,i2, with i ± , ...,ijv £ {t\ 1} form an 

orthonormal basis for Hilbert space Tic- Then, one can easily obtain 
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where, we have used the facts that cr z \ f) = | T))°z| I) — — I I), &x\ T) = I i),Cx| I) — I T) an d 
a y\ T) = A I)) a y\ I) = ~A T)- The equation (3-39) implies that the action of on the basis 
vectors is equivalent to the action of the operator 2P^ — I N , i.e., we have 

H tj = 2P tj - I N , (3-40) 

where, Pij denotes the permutation operator which permutes i-th and j-th sites and In is 
N x N identity matrix, where N is the number of vertices (N := |V|). Now, let |/) denotes the 

vector state which its all components are | except for /, i.e., \l) — | f ... f j | ... |). Then, 

i 

we have 

1 _ _ Nk 

E ^,-10 = j( E p « + 2 E W> = (-2^-«m)IO+ E li>> 

which implies that 

E ^■ = ( /tm( ^~ 2) / + ^)- (3-41) 

(i,j)eRm 

Then, by using (3-40) and (3-41), the hamiltonian in (3-36) can be written in terms of the 
adjacency matrices Ai, i — 0, 1, d as follows 

d d _ 4 

H=J2 J m E ( 2P u - J*) = 2 E J ™^4 m + — r— E ^«™ J - I 3 " 42 ) 

m=0 (i,j)£R m m=0 m=0 

As it has been shown in [27], many known Hamiltonians suitable for PST are basically as- 
sociated with permutations and can thus be obtained within the present unifying theoretical 
framework. For the purpose of the perfect transfer of state, we consider distance-regular 
graphs with n d = \T<i{o)\ = 1, i.e., the last stratum of the graph contains only one site. Then, 
we impose the constraints that the amplitudes {(f)i\e~ lHt \(j) ) be zero for all % — 0, 1, ...,d — 1 
and (0(i|e~ im |0 o ) — e ^ \ where 6 is an arbitrary phase. Therefore, these amplitudes must 
be evaluated. To do so, we use the stratification and spectral distribution associated with 
distance-regular graphs to write 

{(f)i\e- im \(t) Q ) = e - iI ^ li Sl=o JmKm (0 i |e- 2i *^o J ^-|0 o ) = 
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Let the spectral distribution of the graph is fi(x) = Y, k =olk$( x — %k)- Then, (<f>i\e lHt \4>o) — 
implies that 

£ j k P t (x k )e- 2U ^ -W**) = 0, i = 0, 1, d - 1 
Denoting e~ 2rf E m =o J mPm{x k ) ^ above constraints are rewritten as follows 



p i( x k)Vk^k = 0, i = 0,l,...,d- 1, 



fc=0 



(3-43) 



^ P d (x k )r] k -f k = e ld . 

k=0 

As it was discussed previously, Pi(x k ) are entries of the matrix P (P^ = Pi(x k )) which is 
invertible, i.e., the Eq.(3-43) can be written as 



/ \ 




f ° ) 


Volo 













K Vdld ) 




\ 6 J 


rjk'fk for k = 0, 1, ... 


, d are 



(3-44) 



column of the matrix (P*) 1 = ^Q* multiplied with the phase e ld , i.e., for the purpose of PST, 
the following equations must be satisfied 



Vklk = lke- 2Uo ^ JMxk) = —(Q t ) kd , for fc = 0,l,..,d. 



(3-45) 



In the following, we investigate PST between antipodes of some distance-regular networks such 
as cycle networks with even number of nodes and d dimensional hypercube networks. 



3.2 Examples 
1. Cycle graph C 2m 

A well known example of distance-regular networks, is the cycle graph with N vertices denoted 
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by Cat (see Fig. 1 for even N = 2m). For the purpose of perfect transfer of state, we consider 
the cycle graph with even number of vertices, since as it can be seen from Fig. 1, in this case 
the last stratum contains a single state corresponding to the m-th vertex. From Figure 1 it 
can be also seen that, for even number of vertices N = 2m, the adjacency matrices are given 
by 



A = I 2m , Ai^S' + S'', i = \,2,...,m-\, A m = S m , 
where, S is the N x N circulant matrix with period N ( S N = I N ) defined as follows 

^ 1 ... ^ 



(3-46) 



S 



... 10 

... 1 

1 ... 



By using (3-46), one can obtain the following recursion relations for C2m 



(3-47) 



A\Ai = Ai-i + A i+1 , % = 0, 1, ...,m - 1; A\A m — A m _i 



(3-48) 



(the graph C<i m consists of m + 1 strata). By comparing (3-48) with three term recursion 
relations (2-12), we obtain the intersection arrays for C 2m as 



{b , 6 m _i; Ci, c m } — {2, 1, 1, 1; 1, 1, 2}. 



(3-49) 



Then, by using (2-23), the QD parameters are given by 



ai = 0, % = 0, 1, m; uoi = u m = 2, Ui = 1, i = 2, m — 1, 



(3-50) 



By using the recursion relations (2-22), one can show that 



Q {x) = P (x) = 1, Qi(x) = Pi(x) = 2Ti(x/2), i = 1, m-1, Q m (x) = 2P m {x) = 2T m {x/2) 

(3-51) 
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where Tj's are Chebyshev polynomials of the first kind. 

Then, the eigenvalues of the adjacency matrix A = A x (roots of Q m+ \(x) = 2T m+1 (x/2)) 
are given by 

Xi = uj 1 + u)~ l = 2 cos(27ri/iV), % — 0, 1, m 
with uj := e 2m l N . Also, one can show that 7i's (degeneracies of eigenvalues Xi) are given by 



7o = 7m = 1 /2m, 7i = 1/m, i = 1, 2, m - 1. 



(3-52) 



Now, as regards the Eq. (2-5), the matrix P* associated with cycle graph C 2m reads as 



1 1 

2 2 cos(27r/iV) 



1 1 

2cos(2(m- l)n/N) 2u r 



\ 



2 2 cos(2(m - 1)tt/JV) • • • 2 cos((m - ifn/N) 2oj m ( m -^ 
1 w m ••• ••• u m2 



(3-53) 



One can see that (P*) 2 = NI, so the inverse of P* is given by (P*) 1 = ^P*. Therefore, by 



using (3-44) and (3-52), we obtain 



= e -itEi:o 2 -Wos(2WAO) = (-l)V*, i = 0, 1, ...,m 



(3-54) 



For instance, for N = 4, we obtain 



-it (Jo+2Ji+J 2 ) _ e «e 



??i = e -it (Jo-J2) = _ e ie^ 



^ 2 _ e -it {Jo-2.J 1 + .J 2 ) _ 



(3-55) 



which gives us the following equations 



-t{J + 2 Ji + J 2 ) = + 2Ztt, 



-*(Jb - J 2 ) = + (2/' + l)vr, 
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-f (Jo - 2 Ji + J 2 ) = + 2/%. (3-56) 
For I — I' — I" — 0, one can obtain 

= Jl = - j2 = i- (3 - 57) 

whereas by choosing I — I' — 0, 1" — 1, the solution to (3-56) is given by 

Jo = -^, Ji = ^, J 2 = 0. (3-58) 

In the first case, the time to at which the state |0o) = |0) = 1 1000) is perfectly transferd to the 
vertex |0 2 ) = |2) = |0010) is given by 

/ 26 + 71 71 ri^ 

to = -^r = 4^' (3 - 59) 

whereas in the latter case t is given by 

9 + 7T 7T , , 

<» =-—«:■ (3 - 60) 

2. Hypercube network 

The hypercube of dimension d (known also as binary Hamming scheme H(d,2)) is a network 
with N = 2 d nodes, each of which can be labeled by an d-bit binary string. Two nodes on 
the hypercube described by bitstrings x and y are are connected by an edge if \x — y\ = 1, 
where \x\ is the Hamming weight of x. In other words, if x and y differ by only a single bit 
flip, then the two corresponding nodes on the network are connected. Thus, each of 2 d nodes 
on the hypercube has degree d. For the hypercube network with dimension d we have d + 1 
strata with 

The intersection numbers are given by 

bi = d-i, 0<i<d-l; q = i, 1 < i < d. (3-62) 
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Furthermore, the adjacency matrices of this network are given by 



Ai = V cr x <g) g x ... <g) o x ®I 2 <8> ■■■ <8> h, i — 0,l,...,n, 
z — ' » .. .. ' 



(3-63) 



perm. 



where, the summation is taken over all possible nontrivial permutations. In fact, the underlying 
network is the cartesian product of d-tuples of complete network K 2 . Also it can be shown 
that, the idempotents {E , E ly E d } are symmetric product of d-tuples of corresponding 
idempotents of complete network K 2 . That is, we have 



Ei = E - ® E -— ® E - ® E + ® ® E +i i = 0,l,...,tZ, 



(3-64) 



perm. 



d—i 



where 



E ± = -(I±a x ). 



(3-65) 



It has been shown that the eigenmatrices P and Q for the Hamming scheme H(d, 2) are the 
same, i.e., this scheme is self dual [29]. Also, Delsarte [30] showed that the entries of the 
eigenmatrix P = Q for the Hamming scheme H(d, 2) can be found using the Krawtchouk 
polynomials as follows 

Pu = Qu = K l (i), (3-66) 



where Ki(x) are the Krawtchouk polynomials defined as 



Ki(x) = £ 



i=0 



( 



X 



V 1 J 



d — x 



l-i 



:-ir 



(3-67) 



Therefore, we have ({P 



t\-i\ — i 



The eigenvalues xi of the adjacency matrix A = A 1 and corresponding degeneracies 7; are 
given by 

xi — 21 — d] 
d\ 



H = 



2 d l\{d-l)\ 



, Z = 0,l,..,d. 



(3-68) 
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By using (3-66), we have 

Vl = e -^E d m =o J ^ K ^ l \ Z = 0,l,...,d. (3-69) 

Now, in order to evaluate the time t at which PST takes place, the following equations must 
be satisfied 

e i8 e iP 



vm = -^Qdi = 2d K i(d), v 1 = 0,1, d, 



which are equivalent to 

d\ 



e -2it J2 d m ^ J m K m (i) = e * Kl (Q ) y 1 = 0,1, ...,d. (3-70) 



l\(d-l)\ 

For instance, in the case of d = 3 (see Fig. 2), we must solve the following equations 



e -2it (J +3J 1 +3J 2 +J 3 ) _ e i0 



e -2it (J +Jl-J2-J3) _ _ e id 



e -2it (Jo-Ji-J2+J3) — g* 1 



e -2ito(Jo-3Ji+3J2-J 3 ) _ _ e i8 (3-71) 

By solving Eqs. (3-71) one can obtain the following solution 

J ° = J ' = ^ = £ -W, = 0; (3-72, 

that is the time t at which PST takes place is given by 

In the appendix we consider PST over antipodes of some important finite distance-regular 
networks. 
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By using spectral analysis techniques and algebraic combinatoric structures of distance-regular 
graphs (as spin networks) such as stratification introduced in [4, 5] and Bose-Mesner algebra, a 
method for finding a set of coupling constants in some particular spin Hamiltonians associated 
with spin networks of distance-regular type was given so that perfect state transfer between 
antipodes of the networks can be achieved. As examples, the cycle networks with even number 
of vertices and <i-dimensional hypercube networks were considered. 
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Appendix 

In this appendix we consider some important finite distance-regular networks such that their 
last stratum contains only one node. Then by using the prescription of section 3, we investigate 
PST over antipodes of these networks. 
1. Icosahedron [31] 
Intersection array: 

{b , 6i, b 2 ; ci, c 2 , c 3 } = {5, 2, 1; 1, 2, 5}. 
Size of strata and QD parameters: 

Kq — 1, k = Ki — 5, k 2 = 5, /t 3 = 1, 

ao = 0, cci = «2 = 2, 0:3 = 0; a;i = 5, u; 2 = 4, CU3 = 5. 
Polynomials Pj(x): 

P = l, P 1 (x)=x, P 2 (x) = ^(x 2 -2x-5), P 3 (x) = -^(x 3 -4x 2 -5x + 10). 
Stieltjes function: 

x 3 _ 43,2 _ 5a , + 1Q 



GJx) 



M ; x A - 4x 3 - 10x 2 + 20a; + 25' 
Spectral distribution ( = X)f=o7J^( a; ~ x i))'- 

fi(x) = -^{5S(x + 1) + 5(x - 5) + 36(x - Vh) + 35(x + v 7 ^)}. 

Now, one can obtain the matrix P 1 and its inverse. Then by solving the equations (3-45), the 
solution is obtained as follows 

j 69 + 7n (5-3V5> (5 + 3^)71 5tt 

Jo - _ ^T' 1_ 60^' J2 " 60^' J3 -12^- (A "° 

Then, the time t at which PST takes place is given by 
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2. Desargues [31] 
Intersection array: 

{bo, h, b 2 , b 3 , b 4 ; c 1: c 2 , c 3 , c 4 , c 5 } = {3, 2, 2, 1, 1; 1, 1, 2, 2, 3}. 

Size of strata and QD parameters: 

k q = 1, k = Ki — 3, /t 2 = 6, k 3 = 6, k 4 = 3, k 5 = 1 

ctj = 0, i = 0, 1, 5; cji = 3, cj 2 = 2, u> 3 = 4, W4 = 2, CJ5 = 3. 
Polynomials Pi(x): 

P = 1, Pi (a;) = P 2 (x) = a; 2 -3, P 3 (x) = l-(x 3 -5x), P 4 (x) = hx 4 -9x 2 +12), P 5 (x) = -Ua; 5 -lla; 3 +22a;). 
Stieltjes function: 

^ . . x 5 - 11a; 3 + 22a; 

G„ a; = 



x 6 - 14a; 4 + 49a; 2 -36' 
Spectral distribution: 

^(x) = ^r{55(x + 1) + 55(a; - 1) + 45(x + 2) + 45(a; - 2) + 5(x + 3) + 5(x - 3)}. 

The solution to Eq.(3-45) is given by: 

30$ + 517T 7T 47T 7T 7T 

60t 10*o 15t 15*o 4* 

3. Dodecahedron [31] 
Intersection array: 

{6 , 61, 62, &3, ci, c 2 , c 3 , c 4 , c 5 } = {3, 2, 1,1,1; 1, 1, 1, 2, 3}. 
Size of strata and QD parameters: 

Ko — 1, k = Ki = 3, K2 = 6, k 3 = 6, K4 = 3, K5 = 1 
«o = cti = 0, a 2 = a 3 = 1, 0:4 = 0:5 = 0; u;i =3, cj 2 = 2, u; 3 = 1, = 2, cj 5 = 3. 
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Polynomials Pi(x): 

P = l, P 1 (x)=x, P 2 (x)=x 2 -3, P 3 (x) = x 3 -5x-x 2 + 3, P 4 (x) = \{x^ -hx 2 -2x 3 + 

P 5 (x) = ^{x 5 - 7x 3 - 2x 4 + 10x 2 + 10x - 6). 

Stieltjes function: 

x 5 - 7x 3 - 2x 4 + I0x 2 + lOx - 6 



GJx) 



M 1 - 10x 4 - 2x 5 + I6x 3 + 25x 2 - 30x 
Spectral distribution: 

fi(x) = ~ {4<S(a;) + 55{x - 1) + 45{x + 2) + 5(x - 3) + 35{x - Vh) + 35{x + V5)}. 

The solution to Eq.(3-45) is given by: 

9 + 2% (2 + 3V^)n _ 17tt tt _ (2 - 3^)^ 

J °" 2tT' Jl " 60^ ' 2 ""60^' 3 ~60V 4 " 60^ ' 

4. Taylor (P( 13)) [32] 
Intersection array: 

{b , bi, b 2 ; Ci, c 2 , c 3 } = {13, 6, 1; 1, 6, 13}. Size of strata and QD parameters: 

K — 1, = = 13, /t 2 = 13, k 3 = 1, 

«o = 0, cti = «2 = 6, «3 = 0; a;i = 13, u>2 = 36, u 3 = 13. 
Polynomials Pj(x): 

P = l, P 1 (x)=x, P 2 (x) = -{x 2 -6x- 13), P 3 (x) = — (x 3 -12a; 2 -13a; + 78). 

6 78 

Stieltjes function: 

x 3 - 12x 2 - 13x + 78 



x 4 - 12x 3 - 26x 2 + 156x + 169 
Spectral distribution: 

^x) = ^{I35(x + 1) + 5(x - 13) + 75{x - Vl3) + 75{x + ^13)}. 
28 
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The solution to Eq.(3-45) is given by: 

140 + 15tt (13 - 7VT5)tt (13 + 7^)^ _ 13tt 

Jn — , Ja — , Jo — , J'i — . 

28t 364t 364t 28t 

5. Taylor(GQ(2,2)) [32] 
Intersection array: 

{b , b 1 , b 2 ; ci, c 2 , c 3 } = {15, 8, 1; 1, 8, 15}. 
Size of strata and QD parameters: 

Kq — 1, /c = Ki — 15, « 2 = 15, k 3 = 1, 

a = 0, ai = a 2 = 6, a 3 = 0; lo\ — 15, w 2 = 64, cj 3 = 15. 
Polynomials Pi(x): 

P = l, Pi(x)=x, P 2 {x) = \{x 2 -Qx-lh), P 3 (x) = -^-(x 3 -12x 2 -43x + 90). 
Stieltjes function: 

- I2x 2 - 43x + 90 



a Ax) = 



M ; x 4 - 12x 3 - 58x 2 + 180a; + 225' 
Spectral distribution: 

^(x) = — {155{x + 1) + 10<J(x - 3) + 65{x + 5) + <J(x - 15)}. 
32 

The solution to Eq.(3-45) is given by: 

160 + 15tt tt 3tt _ 13tt 

32tg 32tg 32tg 32tg 

6. Taylor(T(6)) [32] 
Intersection array: 

{b , b u b 2 ; ci, c 2 , c 3 } = {15, 6, 1; 1, 6, 15}. 
Size of strata and QD parameters: 

k — 1, /c = Ki = 15, k 2 = 15, k 3 = 1, 
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oo = 0, 0:1 = 0:2 = 8, 03 = 0; lo\ — 15, 002 = 36, 0U3 = 15. 
Polynomials Pi(x): 

P = 1, P^x) = x, P 2 (x) = \{x 2 -8x- 15), P 3 (x) = ^-(x 3 - 16x 2 + 13x + 120). 

o 90 

Stieltjes function: 

„ , , x 3 - 16x 2 + 13x + 120 

GJx) 



M 1 x 4 -16x 3 -2a; 2 + 240x + 225' 
Spectral distribution: 

n(x) = —{I55(x + 1) + 105(x + 3) + 65(x - 5) + 5(x - 15)}. 

The solution to Eq.(3-45) is given by: 

_16#+15tT 37T 7T 137T 

32tg 32tg 32tg 32tg 

7. Wells [33] 
Intersection array: 

{6 , 61, 62, ^3; ci, c 2 , c 3 , c 4 } = {5, 4, 1, 1; 1,1,4, 5}. 
Size of strata and QD parameters: 



Kq — 1, K = Ki = 5, K 2 = 20, K3 = 5, k 4 = 1, 



Oo = 01 = 0, 02 = 3, 03 = 04 = 0; iO\ = 5, 002 = ^3 = 4, co>4 = 5. 
Polynomials Pi(x): 

I 1 

P = 1, Pi(x) = x, P 2 (z) = x 2 -5, P 3 (x) = -(x 3 -9x-3a; 2 +15), P A (x) = —(x A -l?>x 2 -?>x 3 +lhx+2Q). 

Stieltjes function: 

_ a: 4 - 13x 2 - 3x 3 + 15rr + 20 

~ - 18a; 3 - 3x 4 + 30x 2 + 65x - 75 ' 

Spectral distribution: 

^(x) = -Ul0£(a; - 1) + 55(x + 3) + <J(a; - 5) + 85(x + V$) + 85(x - v 7 ^)}. 
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The solution to Eq.(3-45) is given by: 

169 + 23tt _ (5 - 8V5)n 3tt _ (5 + 8^)^ 9tt 
■A) — TT^TL j Jl — 777777 ) ^2 — ~ or> , , J3 — 777777 , Ja 



32t ' 160t ' 32t ' 160t 32t 

8. Hadamard network [34] 
Intersection array: 

{6 , 61, 62, ^3; ci, c 2 , c 3 , c 4 } = {8, 7, 4, 1; 1, 4, 7, 8}. 
Size of strata and QD parameters: 

Kq — 1, k = Ki = 8, « 2 = 14, k 3 = 8, k 4 = 1, 

ccj = 0, % = 0, 1, 4; 0*1 = 8, cj 2 = 28, u; 3 = 28, w 4 = 8. 
Polynomials P{(x): 

P = l, P!(a;)=x, P 2 (a;) = i(a; 2 -8), P 3 (z) = ^(x 3 - 36a:), P 4 (rr) = ^(a; 4 - 64a; 2 + 224). 
Stieltjes function: 

a; 4 - 64a; 2 + 224 
~ x 5 - 72a; 3 + 512a;' 

Spectral distribution: 

^(x) = ^{14,5(3;) + S(x - 8) + + 8) + 85(z - 2^2) + S5(x + 2^2)}. 

The solution to Eq.(3-45) is given by: 

16# + 19tt (1 + 2-v/2)tt 3tt (1 - 2v^)tt _ 13tt 

Jq — , J\ — , t7 2 — , i7 3 — , J 4 — . 

32to 32to 32to 32to 32to 

9. Taylor(P(17)) [32] 
Intersection array: 

{b , bi, b 2 ; Ci, c 2 , c 3 } = {17, 8, 1; 1, 8, 17}. Size of strata and QD parameters: 



K — 1, K = K X = 17, K 2 = 17, K 3 = 1, 
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ao = 0, 0:1 = 0:2 = 8, 03 = 0; lo\ — 17, uj 2 = 64, uj 3 = 17. 
Polynomials Pi(x): 

P = 1, P 1 (x) = x, P 2 (x) = \{x 2 -Bx- 17), P 3 (x) = -^-(x 3 - 16a; 2 - 17 x + 136). 

o 136 

Stieltjes function: 

. * x 3 - I6x 2 - 17x + 136 

GJx) 



" v ' x 4 - 16a; 3 - 34a; 2 + 272a; + 289' 
Spectral distribution: 

fji(x) = ^{l7S(x + 1) + - 17) + 95{x - Vl7) + 95(x + Vl7)}. 

The solution to Eq.(3-45) is given by: 

18# + 19tt (17 - 9-s/l7)7r (17 + 9^)^ _ 17tt 

36t 612t ' 612t 36t 

10. Hadamard network [32] 
Intersection array: 

{6 , 61, b 2 , h; a, c 2 , c 3 , c 4 } = {12, 11, 6, 1; 1, 6, 11, 12}. 

Size of strata and QD parameters: 

Kq — 1, «; = Ki — 12, k 2 = 22, k 3 = 12, /t 4 = 1, 

a, = 0, i = 0, 1, 4; c^! = 12, oj 2 = u 3 = 66, cj 4 = 12. 
Polynomials Pj(x): 

P = l, P 1 (x)=x, P 2 (x) = ^(x 2 -12), P 3 ( x ) = l-( x 3 -78x), P 4 (x) = ^(x 4 - 144x 2 + 792). 
Stieltjes function: 
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Spectral distribution: 

H(x) = ^-{225{x) + 5(x + 12) + 5(x - 12) + 125(x - 2V3) + 125(x + 2^)}. 
48 

The solution to Eq.(3-45) is given by: 

240 + 277T (l-2-s/3)7T 7T (l + 2 V / 3> _ 21tt 

J °~ 48^"' 1_ 48^ ' 2 ~~16V 3 " 48^ ' 4 ~48V 
11. Taylor (SRG (25, 12)) [32] 
Intersection array: 

{b , 6i, 6 2 ; ci, c 2 , c 3 } = {25, 12, 1; 1, 12, 25}. 
Size of strata and QD parameters: 

Kq — 1, = = 25, /t 2 = 25, /t 3 = 1, 

a = 0, cti = a 2 = 12, a 3 = 0; lo\ — 25, cc> 2 = 144, u 3 = 25. 
Polynomials Pi(x): 

P = l, P 1 (x)=x, P 2 ( x ) = -^(x 2 -l2x-25), P 3 ( x ) = ^(x 3 -2Ax 2 - 25x + 300). 
Stieltjes function: 

x 3 - 24a; 2 - 25x + 300 



GJx) 



M ' a; 4 - 24x 3 - 50x 2 + 600x + 625' 
Spectral distribution: 

fi(x) = — {255{x + 1) + 135(x - 5) + 135{x + 5) + 5{x - 25)}. 
52 

The solution to Eq.(3-45) is given by: 

_ _26# + 27tt 2tt 9tt _ 25tt 

°~ 52t^' 1 ~6E^ , 2 ~~130V 3 ~52V 

12. Gosset, Tayl(Schlafli) [32] 

Intersection array: 

{b , 6i, 6 2 ; ci, c 2 , c 3 } = {27, 10, 1; 1, 10, 27}. 
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Size of strata and QD parameters: 

Kq — 1, K = Ki — 27, K 2 = 27, K 3 = 1, 

a = 0, 0:1 = 0:2 = 16, o 3 = 0; LO\ — 27, uj 2 = 100, u 3 = 27. 
Polynomials Pi(x): 

P = l, Pi{x)=x, P 2 {x) = — {x 2 - 16a: -27), P 3 (x) = -^-(x 3 - 32a; 2 + 129x + 432) 

X0 2 (0 

Stieltjes function: 

_ x 3 - 32a; 2 + 129a; + 432 
^ ~ a; 4 - 32a; 3 + 102a; 2 + 864a; + 729' 

Spectral distribution: 

(i(x) = ^{278(x + 1) + 215{x + 3) + 75{x - 9) + 8(x - 27)}. 
00 



The solution to Eq.(3-45) is given by: 



140 + 117T 5n n 5n 

Jo = 7^1 > ^1 = — ' ^2 = j J3 



28t ' 84t 42t ' J Ut 

13. Taylor (Co-Schlafli)) [32] 
Intersection array: 

{60, 61, h; ci, c 2 , c 3 } = {27, 16, 1; 1, 16, 27}. 
Size of strata and QD parameters: 

Kq — 1, /€ = /?! = 27, /t 2 = 27, k 3 = 1, 

oo = 0, oi = 0:2 = 10, 03 = 0; u>i — 27, oj 2 = 256, CJ3 = 27. 
Polynomials Pi(x): 

P = l, P 1 (x)=x, P 2 (x) = -^(x 2 -l0x-27), P 3 (x) = ^(^ 3 - 20a; 2 -183a; + 270) 
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Stieltjes function: 

_ x 3 - 20x 2 - 183x + 270 

~ - 20x 3 - 210x 2 + 540a; + 729' 

Spectral distribution: 

fi( x ) = — {275{x + 1) + 216 (x - 3) + 75(x + 9) + <J(x - 27)}. 
56 

The solution to Eq.(3-45) is given by: 

140 + 117T 7T 5tt 5tt 

^0 — 1 Jl — 7777 ; ^2 — ~ 7777 i ^3 — 



28t 42t 84t 14t 

14. Taylor( 1 Si?G'(29, 14)) [32] 
Intersection array: 

{6 , 61, b 2 ; ci, ca, c 3 } = {29, 14, 1; 1, 14, 29}. 
Size of strata and QD parameters: 

kq — 1, re = «i = 29, re 2 = 29, re 3 = 1, 

«o = 0, «i = «2 = 14, a 3 = 0; c^i = 29, c^2 = 196, 003 = 29. 
Polynomials Pi(x): 

P = l, Pi(x)=x, P 2 (x) = -^(x 2 -Ux-29), P 3 (x) = ^(x 3 - 28x 2 - 29x + 406). 
Stieltjes function: 

x 3 - 28x 2 - 29x + 406 
~ x 4 _ 28x 3 - 58x 2 + 812x + 841 ' 

Spectral distribution: 

fi(x) = -^l 29 ^^ + 1) + S(x - 29) + 156(x - V29) + 155(x + v 7 ^)}. 

The solution to Eq.(3-45) is given by: 

300 + 31tt (29 - 15V29)tt (29 + 15^29)^ _ 29tt 

°~ 60t^' 1_ 1740^ ' 2 ~ 1740^ ' 3 ~60V 
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15. Doubled odd(4) [32] 
Intersection array: 

{&o,&i,&2, h,h,h, b 6 ; ci, c 2 ,c 3 , c 4 ,c 5 , c 6 ,c 7 } = {4,3,3,2,2, 1, 1; 1, 1,2,2,3,3,4}. 

Size of strata and QD parameters: 

Kq — 1, K = Ki = 4, /t 2 = 12, k 3 = 18, «4 = 18, k$ = 12, k 6 = 4, k 7 = 1, 

ttj = 0, i = 0, 1, 7; a;i = 4, u; 2 = 3, cj 3 = 6, u; 4 = 4, cj 5 = 6, cj 6 = 3, w 7 = 4. 
Polynomials Pi(x): 

P = l, Px(x)=a;, P 2 (a;) = a; 2 - 4, P 3 (x) = -(a; 3 - 7x), P 4 (a;) = ~{x 4 - 13a; 2 + 24), 

2 4 

P 5 (x) = -'-(a; 5 -17a; 3 +52a;), P 6 (x) = -Ua; 6 -23a; 4 +130a; 2 -144), P 7 (x) = -^-(a; 7 -26a; 5 +181a; 3 -300a;). 
Stieltjes function: 

- 26a; 5 + 181a; 3 - 300a; 
~ x s _ 30x 6 + 273a; 4 - 820a; 2 + 576' 

Spectral distribution: 

fi(x) = ^{U[5(x-l) + 5(x+l)+5(x-2)+5(x+2)]+6[5(x-3)+5(x+3)]+5(x-A)+5(x+4:)}. 

The solution to Eq.(3-45) is given by: 

706* + 151tt 56tt 47T 71 t /t tt 

^° = 777w ' ^ = °' ^ 2 = ~o/ic;+ ' ^ 3 = °' ^ 4 = mw ' ^ 5 = _ ^T' ^ 6 = ^ 7 = XT' 

140io 245to lOoio 70to oot 4t 

16. J(8,4) [32] 
Intersection array: 

{b , 6i, 62, & 3 ; ci, c 2 , c 3 , c 4 } = {16, 9, 4, 1; 1, 4, 9, 16}. 
Size of strata and QD parameters: 

Kq — 1, = = 16, k 2 = 36, /t 3 = 16, /t 4 = 1, 



Perfect state transfer 34 

«o = 0, Oil = 6, «2 = 8, 0:3 = 6, 0:4 = 0; uj\ — 16, 002 = 36, UJ3 = 36, 0J4 = 16. 
Polynomials Pi(x): 

P = l, P 1 (x)=x, P 2 (x) = ^(x 2 - 6x - 16), P 3 (^) = ^(a: 3 -14a; 2 -4a;+128), 
P 4 (x) = -^(x 4 - 20a; 3 + Ux 2 + 368x - 192). 

Stieltjes function: 

„ , . x 4 - 20x 3 + 44x 2 + 368a; - 192 

GJx) 



M 1 - 20x 4 + 28x 3 + 592a: 2 -128a; - 2048' 
Spectral distribution: 

H(x) = ^{S(x - 16) + 75(x - 8) + U5(x + 4) + 28<5(a; + 2) + 205(a; - 2)}. 

The solution to Eq.(3-45) is given by: 

706 + 199vr 7t 13tt tt 17tt 

^0 — 7777: ) <A — or . ) ^2 — 777777"! ^3 — — . , , J 4 — 



140t 35t 210t 14*o 140t 
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Figure Captions 

Figure 1: Denotes the cycle network C 2m , where the m+1 vertical dashed lines show the 
m + 1 strata of the network. 

Figure 2: Shows the cube or Hamming scheme if (3, 2) with vertex set V = {(ijk) : 
i, j, k = 0, 1}, where the vertical dashed lines denote the four strata of the cube. 



